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Abstract — The growth rate of the weight distribution of irreg- 
ular doubly-generalized LDPC (D-GLDPC) codes is developed 
and in the process, a new efficient numerical technique for 
its evaluation is presented. The solution involves simultaneous 
solution of a 4 x 4 system of polynomial equations. This represents 
the first efficient numerical technique for exact evaluation of the 
growth rate, even for LDPC codes. The technique is applied to 
two example D-GLDPC code ensembles. 



I. Introduction 

Recently, the design and analysis of coding schemes repre- 
senting generalizations of Gallager's low-density parity-check 
(LDPC) codes [ 1 ] has gained increasing attention. This interest 
is motivated above all by the potential capability of these 
coding schemes to offer a better compromise between waterfall 
and error floor performance than is currently offered by state- 
of-the-art LDPC codes. 

In the Tanner graph of an LDPC code, any degree-q variable 
node (VN) may be interpreted as a length-q repetition code, 
i.e., as a (q, 1) linear block code. Similarly, any degree-s check 
node (CN) may be interpreted as a length-s single parity- 
check (SPC) code, i.e., as a (s, s — 1) linear block code. The 
first proposal of a class of linear block codes generalizing 
LDPC codes may be found in [2|, where it was suggested 
to replace each CN of a regular LDPC code with a generic 
linear block code, to enhance the overall minimum distance. 
The corresponding coding scheme is known as a regular 
generalized LDPC (GLDPC) code, or Tanner code, and a CN 
that is not a SPC code as a generalized CN. More recently, 
irregular GLDPC codes were considered (see for instance 0). 
For such codes, the VNs exhibit different degrees and the CN 
set is composed of a mixture of different linear block codes. 

A further generalization step is represented by doubly- 
generalized LDPC (D-GLDPC) codes |g). In a D-GLDPC 
code, not only the CNs but also the VNs may be represented by 
generic linear block codes. The VNs which are not repetition 
codes are called generalized VNs. The main motivation for 
introducing generalized VNs is to overcome some problems 
connected with the use of generalized CNs, such as an overall 
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code rate loss which makes GLDPC codes interesting mainly 
for low code rate applications, and a loss in terms of decoding 
threshold (for a discussion on drawbacks of generalized CNs 
and on beneficial effects of generalized VNs we refer to [5] 
and (6), respectively). 

A useful tool for analysis and design of LDPC codes and 
their generalizations is represented by the growth rate of the 
weight distribution, or equivalently, the asymptotic weight 
enumerating function (WEF). The growth rate of the weight 
distribution was introduced in HI to show that the minimum 
distance of a randomly generated regular LDPC code with a 
VN degree of at least three is a linear function of the codeword 
length with high probability. The same approach was taken in 
and ID to obtain related results on the minimum distance 
of subclasses of Tanner codes. 

The growth rate of the weight distribution has been subse- 
quently investigated for unstructured ensembles of irregular 
LDPC codes. Works in this area are J9), ED, 03], Ell- 
in particular, in lfl2l a technique for approximate evaluation 
of the growth rate of any (eventually expurgated) irregular 
LDPC ensemble has been developed, based on Hayman's 
formula. Asymptotic weight enumerators of ensembles of 
irregular LDPC codes based on protographs and on multiple 
edge types have been derived in [13] and fl4l . respectively. 
The approach proposed in |[T3"1 has then been extended to 
protograph GLDPC codes and to protograph D-GLDPC codes 
in 1(151 and Ifl6l , respectively. In IfTTl , the authors presented a 
compact formula for the growth rate of general unstructured 
irregular D-GLDPC code ensembles for the specific case of 
small weight codewords. 

In this paper, an analytical expression for the growth rate 
of the weight distribution of a general unstructured irregular 
ensemble of D-GLDPC codes is developed. As opposed to the 
formula developed in |[T7l . the proposed expression holds for 
any codeword weight. The present work also extends to the 
fully-irregular case an expression for the growth rate obtained 
in ifTHl assuming a CN set composed of linear block codes 
all of the same type. In the process of this development, we 
obtain an efficient evaluation tool for computing the growth 
rate exactly. This tool always requires the solution of a (4 x 4) 
polynomial system of equations, regardless of the number 
of VN types and CN types in the D-GLDPC ensemble. 
As shown through numerical examples, the proposed tool 
allows to obtain a precise plot of the growth rate with a low 
computational effort. For the case of irregular LDPC codes, 
a technique for numerical evaluation of the growth rate of 



the weight distribution was given in [12]; in contrast to the 
technique developed in this paper, the method of [12] provided 
an approximate numerical solution for the growth rate; it is 
also more computationally complex than that proposed in the 
present work. 

II. Preliminaries and Notation 

We define a D-GLDPC code ensemble M. n as follows, 
where n denotes the number of VNs. There are n c different 
CN types t £ I c = {1, 2, • • • ,n c }, and n v different VN 
types t £ I v — {1,2,- ■• ,n v }. For each CN type t £ I c , 
we denote by h t , St and r t the CN dimension, length and 
minimum distance, respectively. For each VN type t £ I v , 
we denote by k t , qt and p t the VN dimension, length and 
minimum distance, respectively. For t £ I c , pt denotes the 
fraction of edges connected to CNs of type t. Similarly, for 
t £ I v , At denotes the fraction of edges connected to VNs of 
type t. Note that all of these variables are independent of n. 

The polynomials p(x) and X(x) are defined by p(x) = 
Ete/ C Ptx 8 ^ 1 and A(x) = J2tei v ^x qt ~ 1 - If E denotes the 
number of edges in the Tanner graph, the number of CNs of 
type t £ I c is then given by Ept/st, and the number of VNs 
of type t £ I v is then given by EXt/qt- Denoting as usual 
L p(x) &x and X(x) dx by J p and J A respectively, we 
see that the number of edges in the Tanner graph is given by 
E = n/ J X and the number of CNs is given by m = E J p. 
Therefore, the fraction of CNs of type t £ I c and the fraction 
of VNs of type t £ I v are given by 

it = and 5 t = ~tt W 

s t J p qt jX 

respectively. Also the length of any D-GLDPC codeword in 
the ensemble is given by 

EXt 

qt 
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E 

te/„ 



n x - X t k t 

~ x tk~ 



(2) 



Note that this is a linear function of n. Similarly, the total 
number of parity-check equations for any D-GLDPC code in 
the ensemble is given by M = jj J2tei c pt( ^ fet) . A mem- 
ber of the ensemble M. n then corresponds to a permutation 
on the E edges connecting CNs to VNs. 

The WEF for CN type t £ I c is given by A^(z) = 1 + 
J2u=r t A^u 1 z w . Here > denotes the number of weight- 
it codewords for CNs of type t. The input-output weight 
enumerating function (IO-WEF) for VN type t £ I v is given 
by y) = l + EtU Ev= Pt B { uU u y v . Here > 

denotes the number of weight-w codewords generated by input 
words of weight u, for VNs of type t. Also, is the total 
number of weight-2 codewords for VNs of type t. 

If there exist CNs and VNs with minimum distance equal 
to 2, and define the (positive) parameters 



C = 2 £ ^;V = 2 V- 



St 



t : pt=2 



qt 



(3) 



t : r t =2 

The design rate of any D-GLDPC ensemble is given by 

E teIc Pt(l - Rt) 



R = l 



(4) 



where for t £ I c (resp. t £ I v ) R t is the local code rate of a 
type-* CN (resp. VN). 

Throughout this paper, the notation e = exp(l) denotes 
Napier's number, all the logarithms are assumed to have base 
e and for < x < 1 the notation h(x) = — xlog(x) — (1 — 
x) log(l — x) denotes the binary entropy function. 

III. Growth Rate of the Weight Distribution of 
General Irregular D-GLDPC Code Ensembles 

The growth rate of the weight distribution of the irregular 
D-GLDPC ensemble sequence {.M„} is defined by 



G(a) = lira -logE^,,, [N a 

n — >oo fi 



(5) 



where E^vi„ denotes the expectation operator over the ensem- 
ble M n , and N w denotes the number of codewords of weight 
w of a randomly chosen D-GLDPC code in the ensemble M. n . 
The limit in <(5j assumes the inclusion of only those positive 
integers n for which an £ Z and E^4 n [N an ] is positive. Note 
that the argument of the growth rate function G(a) is equal 
to the ratio of D-GLDPC codeword length to the number of 
VNs; by (0, this captures the behaviour of codewords linear 
in the block length, as in [12] for the LDPC case. 

A D-GLDPC ensemble is said to be asymptotically good 
if and only if a* = inf{a > | G{a) > 0} > 0. 
The parameter a* is called the ensemble relative minimum 
distance. In |[T9l , it was shown that a D-GLDPC ensemble is 
always asymptotically good if there exist no CNs or VNs with 
minimum distance 2 while, if the exist both CNs and VNs with 
minimum distance 2, the ensemble is asymptotically good if 
and only if C ■ V < 1, where C and V are given by (01. 

Note that using (0, we may also define the growth rate with 
respect to the number of D-GLDPC code bits N as follows: 



ff( 7 )= lim -\ogE Mn [N lN } 

N— >-oo iv 



It is straightforward to show that 



H(i) 



(6) 



(7) 



where 



1 \ - Xth 



In this section, we formulate an expression of the growth 
rate for an irregular D-GLDPC ensemble M. n over a wider 
range of a than was considered in lfT7l . |fl9l (where the case 
a — > was analyzed). 

The following theorem constitutes our main result. 

Theorem 1: The growth rate of the weight distribution of 
the irregular D-GLDPC ensemble sequence {.M ra } is given by 



G (a) = ^logB !t) (^Oiyo) ~ a\ogx 



telv 



{f)E7 S lo g ^(zo) + 1 ° S(l - ; f /A) (8) 



where xq, yo, Zq and are the unique positive real solutions 
to the 4x4 system of polynomial equation^ 



JjP 
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f 1 BW(x 07 y ) 



yo s t 
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and 



- c s(')(xo,yo) 

/3 / A ) (1 + y z ) = y Q z a 



(9) 



(10) 



(11) 



(12) 



The theorem is proved in Section [TV] 

IV. Proof of the Main Result 

In this section we prove Theorem Q] The proof uses the 
concepts of assignment and split assignment, defined next. 

Definition 1: An assignment is a subset of the edges of the 
Tanner graph. An assignment is said to have weight k if it 
has k elements. An assignment is said to be check-valid if the 
following condition holds: supposing that each edge of the 
assignment carries a 1 and each of the other edges carries a 
0, each CN recognizes a valid local codeword. 

Definition 2: A split assignment is an assignment, together 
with a subset of the D-GLDPC code bits (called a codeword 
assignment). A split assignment is said to have split weight 
(u, v) if its assignment has weight v and its codeword as- 
signment has u elements. A split assignment is said to be 
check-valid if its assignment is check-valid. A split assignment 
is said to be variable-valid if the following condition holds: 
supposing that each edge of its assignment carries a 1 and 
each of the other edges carries a 0, and supposing that each 
D-GLDPC code bit in the codeword assigment is set to 1 and 
each of the other code bits is set to 0, each VN recognizes a 
local input word and the corresponding valid local codeword. 

For ease of presentation, the proof is broken into two parts. 

A. Number of check-valid assignments of weight 8m 

First we derive an expression, valid asymptotically, for the 
number of check-valid assignments of weight 8m. For each 
t € I c , let etm denote the portion of the total weight 6m 
apportioned to CNs of type t. Then e t > for each t G I c , 
and 2~2tei e * = ^- ^ so denote e = (ej £2 ■ •• e n c )- 

Consider the set of jtm CNs of a particular type t G I c , 
where j t is given by (Q~|). Using generating functions, the 
number of check-valid assignments (over these CNs) of weight 
e t m is given by 

7 t m 



N, 



c.t 



(e t m) = Coeff (A (t) (x) 



,x 



'Note that while {5), )10t and {TT} are not polynomial as set down here, 
each may be made polynomial by multiplying across by an appropriate factor. 



where Coeff [p(x),x c ] denotes the coefficient of x c in the 
polynomial p(x). We next make use of the following result, 
which is a special case of lfl2l Corollary 16]: 

Lemma 1: Let A(x) = 1 + X)u=c A U X V , where 1 < c < d, 
be a polynomial satisfying A c > and A u > for all c < 
u < d. Then 



lim \ log Coeff 

1^00 i 



{A(x))\x^ 



log 



A(z) 



where z is the unique positive real solution to 

A'(z) 



A(z) 



(13) 



(14) 



Applying this lemma by substituting A(x) = A^\a 
jtni and £ = £4/74, we obtain that as m — > 00 

7t«i 



iV c ( ] tm) (e t m) = Coeff 



exp 



jm ^7 t log^ (t) (z 0: t) - et logz ,t)} 



(15) 
(16) 



where, for each t G I c , z o t is the unique positive real solution 
to 



7* 



zo.t = e t 



(17) 



The number of check-valid assignments of weight 6m 
satisfying the constraint e is obtained by multiplying the 
numbers of check-valid assignments of weight e t m over 7 t m 
CNs of type t, for each t G I c , 

= n<rW). 

tela 



Nf\6r 



(18) 



The number N c (8m) of check- valid assignments of weight 
8m is then equal to the sum of 7V c (e) (6m) over all admissible 
vectors e; therefore by ( fl6] >, as m — * 00 



N c (8m) -> ^2 exp{mW(e)} 



where 



W ( e ) = £ (7 t log^ (t) U, t )-etlog 



zo.t 



tela 



(19) 



(20) 



As m — ► 00, the asymptotic expression is dominated by the 
distribution e which maximizes the argument of the exponen- 
tial function^ Therefore as m — > 00 



N c (6m) — y exp {mX} 



where 



X = maxVF(e) 

e 

and the maximization is subject to the constraint 

V(e) = Y,tt = 6 
tela 



(21) 
(22) 

(23) 



together with e f > for each t G I and for every t G I c , 
zo.t is the unique positive real solution to ([T7| >. Note that for 

2 Observe that as m — > 00, exp(mZt) — > exp(mmaxt{Zt}) 



each t G I c , dTTt provides an implicit definition of zo,t as a 
function of e f . 

We solve this optimization problem using Lagrange multi- 
pliers, ignoring for the moment the inequality constraints. At 
the maximum, we must have 

for all t £ I c , where A is the Lagrange multiplier. This yields 

dz , t . — y— i-n./.) _e^_ 

Zo,t 



9e t 



7t 



A(*)(2b )t ) 



log zo„ 



A 



(25) 



The term in square brackets is equal to zero due to (F7\ ; 
therefore this simplifies to logzo.t = —A for all t 6 I c . We 
conclude that all of the {zq,*} are equal, and we may write 



(26) 



we obtain 



zo,t = zo yt e I c . 

Making this substitution in (f2Tb and using 

N c (5m) exp jm ^JT j t \og A^(z Q ) - <51ogz ^ | . 

(27) 

Summing (fTTI i over t € I c and using d23l > and (l26*l l implies 
that the value of zq in (|27| | is the unique positive real solution 
to © (here we have also used the fact that n J p = m / A). 

B. Polynomial-System Solution for the Growth Rate 

Consider the set of <$ t n VNs of a particular type i G J„, 
where 5 t is given by (Q3. Using generating functions, the 
number of variable-valid split assignments (over these VNs) 
of split weight (a t n,(3 t n) is given by 

N$ n) (a t n, (3 t n) = Coeff (flW (x, y)) *"* , 

where Coeff [p(x, y), x c y d ] denotes the coefficient of x c y d in 
the bivariate polynomial p(x, y). We make use of the following 
result, which is a special case of [12, Corollary 16]: 

Lemma 2: Let 

k d 
u—1 v—c 

where k > 1 and 1 < c < d, be a bivariate polynomial 
satisfying £?„ „ > for all 1 < u < k, c < v < d. Then 



lim ilogCoeff (B(x,y)) i ,x^y ei = log , 

L \ 4vo 

where xq and yo are the unique positive real solutions to the 
pair of simultaneous equations 



/ B(x ,yo) 



and 



S(a; ,yo) 



• yo = 



(29) 



(30) 



Applying this lemma by substituting B(x,y) = B^(x,y), 
I — Stn, £ = at/St and 6 = /3 t /<$ t , we obtain that as n — > oo 

(S £ n 



N^ n \a t n,(3 t n) = Coeff 



exp 



{nX t (5t) (a t ,/3 t )} (31) 



where 

X t (l5t) (a t , fa) = S t log B w (x ,t,yo,t)-a t log x ,t-fit log j/o,t 

(32) 

and where xo.t and yo,t are the unique positive real solutions 
to the pair of simultaneous equations 



(j:r 



■(xo,t,Vo,t) 



BW{x ,t,yo,t) 



■ x ,t = a> t 



and 



, -d^-{xo,t,yo,t) 

Ot-^rrr, r ■ Vo,t = Pt 



(33) 



(34) 



B(t)(x ,uyo,t) 

Next, note that the expected number of D-GLDPC code- 
words of weight an in the ensemble A4 n is equal to the 
sum over of the expected numbers of split assignments of 
split weight (an, j3n) which are both check-valid and variable- 
valid, denoted N^ 0n : 

E Mri [N an }=Y / ^M n [NZ 0n }. 



This may then be expressed as 



E E 



Pc 



valid 



{fin) 



a t >o,tei v 0t>o,tei v 

J2 t a t =a 



Y[N$ n \atn,ptn). (35) 



t<=i v 

where (3 — 2~2tei @t- Here P c -vaiid(/9^) denotes the probability 
that a randomly chosen assignment of weight j3n is check- 
valid, and is given by 

Pc-vaiid(/3n) = N c ((3n)/(p n 

Applying (El eqn. (25)], we find that as n — » oo 

K /3nJ \ j3n 
Combining this result with d27l ), we obtain that as n — > oo 

-valid (f3n) — > exp {nY (/?)} 

where 

7pW.,„J,( t )^ h(p J a) 



exp <{ y-^/i ( /3 / A 



^ E^ lo g(^ (t) (-o))-/31og 2 o-^ 



A 



Therefore, as n — > oo 

Eai„ [V Qn ] -f 



E E exp i n 

a t >o,teJ„/3 t) iei"„ I \*e/„ /J 



(36) 



where 



(37) 



tei v 



Note that the sum in Q6j is dominated asymptotically by the 
term which maximizes the argument of the exponential func- 
tion. Thus, denoting the two vectors of independent variables 

by a = (a t ) teIv and j3 = (fi t )tei v , we have 



G(a) = m&xS(ot, P) 

a,/3 



(38) 



Y({3) 



(39) 



where 

tei v 

where fi is given by (l37k and the maximization is subject to 
the constraint 

R(a,/3) = ^2a t = a (40) 



together with a t > and appropriate inequality constraints on 
fit for each t G I v , and J2t a t = a - 

Note that (O provides an implicit definition of zq as a 
function of (3. Similarly, for any t G I v , ( 1331 and (l34l provide 
implicit definitions of XQ, t and yo,t as functions of the two 
variables at and /3t. 

We solve the constrained optimization problem using La- 
grange multipliers, ignoring for the moment the inequality 
constraints. At the maximum, we must have 

dS(a,/3) _ dR{a,/3) 
dat ' dat 
for all t G I v , where /i is the Lagrange multiplier. This yields 



dx -, 
da t 



SB 



'—(xo,t,Vo,t) a t 



S(*)(x ,t,yo,t) 5Co,t 



log x Qi , 



+ 



dyo,t 

da t 



<)B 



dZ-(xo,t,yo,t) fit 



B (t 'Kt,!/o,t) 2/o,t 



The terms in square brackets are zero due to (l33l and (l34l 
respectively; therefore this simplifies to log so, t = — fJ, for all 
t E I v . We conclude that all of the {xo.t} are equal, and we 
may write 

x ,t =x Vt G 4 . (41) 
At the maximum, we must also have 

dS(a,/3) _ dR(a,p) 



9/3t 

for all t E I v . This yields 



dx j 
da t 

dy , t 
dfit 



.IB 



'—{xo,t,yo,t) a t 



B( t} (xo,t,yo,t) xo,t 
-5j-{xo,t,y ,t) fit 



dz 
dfit 



Bi t ){xQ, t ,yo,t) yo,t 



logj/o.t - log 20 
-fij\ 



log 



/A 



E 

se/ c 



7s 



(*o) /? 



fiJX 



. (42) 



The terms in square brackets are zero due to d33l l, 041 ) and 
(O respectively; therefore this simplifies to 

1-fifX' 



zoyo,t 



fiJX 



= 1 Vtel v . 



(43) 



We conclude that all of the {yo,t} are equal, and we may write 

Vo,t = yo Vtelc (44) 

Rearranging j43l we obtain ( TT2b . Also, summing ( f33l > over 
t E I v and using d40t and (|4TT > yields ( fTOb . Similarly, summing 
d34l > over t G I v and using d37j and (0U) yields dTTb . 
Substituting back into ([39]l and using gT}, d44l . dtol l and (l37b 
yields 

= E 5 * iog- 8 ^ (^0,2/0 ) -aloga;o -/31ogy 



JjP 



] ^7 s logA( s )(z )-/31og2 



/ A 



(45) 



where xq, yo, z$ and fi are the unique positive real solutions to 
the 4 x 4 system of equations (0, ([Tol l. ( ITTb and ( fT2l . Finally, 
(fT2l leads to the observation that 



-/31ogz -^logyo 



h(fifX) _ log (1-/3/ A) 



/ A /A 
which, when substituted in ( f45l >. leads to (O. 

V. Examples 

In this section the growth rates of two example D-GLDPC 
ensembles of design rate R = 1/2 are evaluated using the 
polynomial solution of Theorem [TJ We use Hamming (7, 4) 
codes as generalized CNs and SPC codes as generalized VNs. 
Three representations of SPC VNs are considered, namely, 
the cyclic (C), the systematic (S) and the antisystematic (A) 
representations 0. 

Ensemble 1 is characterized by two CN types and two 
VN types. Specifically, we have I c = {1,2}, where 1 G I c 
denotes a (7, 4) Hamming CN type and 2 £ J c denotes a 
length-7 single parity check (SPC) CN type, and = {1, 2}, 
where 1 G I v denotes a repetition-2 VN type and 2 G I v 
denotes a length-7 SPC CN type in cyclic form. Ensemble 
2 is characterized by two CN types and four VN types. 
Specifically, we have I c = {1,2}, where 1 G J c denotes 
a (7,4) Hamming CN type and 2 G I c denotes a SPC-7 
CN type, and I v = {1,2,3,4}, where 1 G I v denotes a 
repetition-2 VN type, 2 G I v denotes a length-7 SPC CN 
type in cyclic form, 3 G I v denotes a length-7 SPC CN 
type in antisystematic form, and 4 G denotes a length-7 
SPC CN type in systematic form. The edge-perspective type 
distributions of the two ensembles are summarized in Table U 

Both Ensemble 1 and Ensemble 2 have been obtained by 
performing a decoding threshold optimization with differential 
evolution (DE) [20|. Ensemble 1 has been obtained by only 

3 The (k X (k + 1)) generator matrix of a SPC code in A form is obtained 
from the generator matrix in S form by complementing each bit in the first 
k columns. Note that a (k X (k + 1)) generator matrix in A form represents 
a SPC code if and only if the code length q = k + 1 is odd. For even k + 1 
we obtain a d mm = 1 code with one codeword of weight 1. 



TABLE I 

Coefficients of X(x) and p(x) for the two example D-GLDPC 

ENSEMBLES. 



presented for two example optimized irregular D-GLDPC code 
ensembles. 



Ensemble 1 


Variable nodes 




Check nodes 






Lrepetition— 2 


Ai 


= 0.055646 l:Hamming(7,4) 


Pi 


= 0.965221 


2:SPC-7 (C) 


A 2 


= 0.944354 2:SPC-7 


P2 


= 0.034779 


Ensemble 2 


Variable nodes 




Check nodes 






l:repetition— 2 


Ai 


= 0.022647 l:Hamming(7,4) 


Pi 


= 0.965221 


2:SPC-7 (C) 


A_> 


= 0.100000 2:SPC-7 


P2 


= 0.034779 


3:SPC-7 (A) 


A 2 


= 0.539920 
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Fig. 1. Growth rates of the two example ensembles described in Section Fvl 
Ensemble 1 is asymptotically bad, while Ensemble 2 is asymptotically good 
with an ensemble relative minimum distance of a* = 2.625 X 10 — 3 . 



imposing the node type and R — 1/2 constraint. In this case 
we have C ■ V = 1.19 > 1, so the ensemble is asymptotically 
bad (a* = 0). Ensemble 2 has been obtained by imposing 
the node type and R = 1/2 constraint, together with the 
constraints C ■ V < 0.5 and A2 > 0.1. Since in this case we 
have C • V = 0.5 < 1, the ensemble is asymptotically good 
(a* > 0). The expected asymptotically bad or good behavior 
of the two ensembles is reflected in the growth rate curves 
shown in Fig. Q] Using a standard numerical solver, it took 
only 5.1 s and 6.7 s to evaluate 100 points on the Ensemble 1 
curve and on the Ensemble 2 curve, respectively. The relative 
minimum distance of Ensemble 2 is a* — 2.625 x 10 -3 . 

VI. Conclusion 

A general expression for the asymptotic growth rate of the 
weight distribution of irregular D-GLDPC ensembles has been 
presented. Evaluation of the expression requires solution of 
a 4 x 4 polynomial system, irrespective of the number of 
VN and CN types in the ensemble. Simulation results were 
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